
setup;

A = prepareA;
f = prepareF;

v = solver(A, f);

% maybe reshape here?..
surf(v);
surf(abs(v-u));
